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We study the Mott phases and the superfluid-insulator transition of two-component ultracold 
bosons on a square optical lattice in the presence of a non-Abelian synthetic gauge field, which 
renders a SU (2) hopping matrix for the bosons. Using a resummed hopping expansion, we calculate 
the excitation spectra in the Mott insulating phases and demonstrate that the superfluid-insulator 
phase boundary displays a non-monotonic dependence on the gauge field strength. We also compute 
the momentum distribution of the bosons in the presence of the non-Abelian field and show that they 
develop peaks at non-zero momenta as the superfiuid-insulator transition point is approached from 
the Mott side. Finally, we study the superfiuid phases near the transition and discuss the induced 
spatial pattern of the superfiuid density due to the presence of the non-Abelian gauge potential. 



I. INTRODUCTION 

Experimental systems involving ultracold bosons in 
optical lattices provide us vifith a unique test bed for 
studying quantum phase transitions [1] |2] and for mim- 
icking strongly correlated condensed matter systems [3]. 
It is well-known that the simplest of such systems with 
bosonic ultracold atoms are well represented by the Bose- 
Hubbard (BH) model for which the superfluid-insulator 
transition has been studied both theoretically [4H9] and 
experimentally pQ [TO]. These ideas have also been ex- 
tended to systems of multi-species and higher spin bosons 
which provide richer phase diagrams [llH13j . Such theo- 
retical works have been supported by recent experiments 
involving studies of both equilibrium phase transitions 
and closed non-equilibrium dynamics in some of these 
systems [l4l[T5] . 

More recently, there have been several theoretical pro- 
posals and concrete experimental realization of synthetic 
Abelian vector potentials for ultracold gases by subject- 
ing them to carefully tuned Raman lasers (for a review 
see Ref. [16]). These lasers induce a space or time de- 
pendent shift in the energy dispersion of the atoms at 
a given momenta and thus emulate the effect of a syn- 
thetic magnetic/electric field [171 - 123] . It has been shown 
that spatially varying Abelian gauge potentials acting on 
a bosonic gas within an optical lattice have a profound 
effect on both the superfluid-insulator transition and the 
superfluid (SF) state to which the transition takes place 
pl'-^^F. In particular, the positions of the precursor Mott 
peaks shift to flnite momenta and the SF density develops 
a spatial ordering pattern which can be detected in stan- 
dard experiments. Similar, albeit more complicated, the- 
oretical proposals have been put forward for realization of 
synthetic SU(2) gauge potentials for neutral atoms 
[2^ . However, the effect of non-Abelian gauge potentials 
on the superfluid-insulator transition of ultracold bosons 
has not been studied so far. 

In this paper we consider a two-species boson sys- 



tem described by a BH model in a square optical lat- 
tice and study the effect of a SU (2) gauge potential on 
its superfluid-insulator transition. Our work involves an 
extension of the resummed hopping expansion in Ref. 
l30l and [31] and is somewhat similar to the strong cou- 
pling expansion technique of Ref. [7] Here these meth- 
ods are applied to the case of two-species bosons in the 
presence of an SU{2) gauge fleld. The main results of 
this work are the following. First, we show that even 
for the simplest non-Abelian gauge fleld, described by a 
constant SU{2) gauge potential, the interplay of inter- 
and intra-species interactions and the presence of the ex- 
ternal gauge fleld leads to qualitative changes in several 
aspects of the superfluid-insulator transition. Second, 
we compute the momentum distribution of the bosons in 
the Mott phase near the superfluid-insulator transition 
and show that the precursor peaks occur at flnite, rather 
than zero, momenta due to the presence of the gauge 
field. We find a sudden change of the peak positions 
when the strength of the gauge field is varied in certain 
parameter regimes. This indicates that a slight change 
in the gauge field strength has strong impact on the dy- 
namical behavior of the system and might be especially 
relevant for its non-equilibrium dynamics. This behav- 
ior is somewhat reminiscent of the QPT of excited states 
discussed in Ref. and of the abrupt sign change of the 
Hall conductivity which has recently been found for a 
system of hardcore bosons |33j . Third, we demonstrate 
that the superfluid-insulator phase boundary displays a 
non-monotonic dependence on the strength of the gauge 
field leading to re-entrant superfluid-insulator transitions 
with the variation of the strength of the gauge fleld for 
a fixed hopping strength. Finally, we construct an effec- 
tive Landau-Ginzburg theory for the superfluid-insulator 
transition and use it to chart out the nature of the SF 
phase into which the transition takes place. We show that 
for a constant non-Abelian gauge fleld, the SF density 
near the transition does not exhibit any spatial ordering. 
This feature is to be contrasted with the case where an 
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Abelian flux (with half flux quanta per lattice plaquette) 
is added over the existing SU{2) potential leading to a 
spatial pattern in the SF density. 

The organization of the rest of the paper is as follows. 
We outline the model describing our system in Section 
nl The Mott insulating (MI) phase is analyzed in Section 
Tn) This is followed by the analysis of the SF phase in 
Section |IV| Finally, we conclude in Sec. [V} 



II. MODEL 

It is well-known that a SU (2) gauge field can be gener- 
ated for a system of two-species ultracold bosons by dis- 
tinguishing between the atoms in two different Zeeman 
levels representing the two 'flavors' of the non-Abelian 
theory [27] ■ In doing so, one substitutes the standard 
hopping process by a laser-assisted tunneling, which may 
depend on the position and the state of the atom. One 
of the crucial features of such a tunneling is its ability to 
flip the state of the atom. Taking into account local in- 
teractions between the atoms, the effective Hamiltonian 
describing the system can be seen to be the same as that 
of a two-species BH model [TTl [T^ |^, with an addi- 
tional non-Abelian vector potential in the hopping term 
providing the additional inter-species coupling. Defining 
the number operators fif = aJiiki and fi^ = bjbi, where 
ki and hi denote the boson annihilation operators of the 
two species, the local part of the Hamiltonian reads 



Hn 



E [-^^tii^t - 1) + ^fi'(n^ - 1) (1) 



[;ab^a-b _ ^a-a 



where U^^ is the strength of interactions between a pair 
of particles with flavors x and y and ii^ denotes the chem- 
ical potential of species x. The Hamiltonian Ho is eas- 
ily diagonalized using the Fock state basis: ]io\n'^,n^) = 
„b jn'^, v}'). It is easy to see that Hq allows for ground- 
state degeneracies which are lifted in the presence of a 
hopping term leading to different types of "magnetic" 
orderings in the Mott state [TT]. 

As we wish to focus on the influence of gauge fields 
on the MI-SF transition, wc shall first consider the pa- 
rameter regime for which Hq has a unique ground-state 
in the local limit. The simplest choice in this regard 
is to assume two independent, identical systems, i.e. 
f^aa ^ jjbh = u , fi'' = fj.^ = fi, and U""^ = 0. With this, 
Hq describes a system which in its ground-state is in both 
flavors occupied by an equal number n of particles, with 
fi/U < n < fi/U + 1 as in the one-component BH model. 
This setup will be further investigated in Sections |IIIB| 
and |IV[ for the MI and SF phase, respectively. 

In another scenario, investigated in Section |III C and 
|IV[ we take into account repulsive interactions between 



the components {U > 0). In this case, a subspace 



degenerate ground state manifold, where h is the parti- 
cle number per site. For < fi/U < If^^/U, one single 
atom occupies each site, such that any linear combina- 
tion ci |1, 0) -t-C2|0, 1) at every site is a ground state of the 
unperturbed Hamiltonian Hq. The hopping lifts this de- 
generacy. As shown in Ref. [121 for small A = U""^ /U , an 
antiferromagnetic ordering is preferred, while for A ~ 1, 
the system chooses a translational-invariant ferromag- 
netic phase. Between these two limits an XY ordering 
with ci = C2 = l/v^ occurs. In the rest of this work, 
we shall focus on the ferromagnetic and the XY phases, 
where the ground state preserves translational symmetry. 
The kinetic part of the Hamiltonian is given by 



(2) 



where = 5<ij>Je ^ 



is a nearest-neighbor 
hopping with a constant strength J, and we have chosen 
h — c—1. The phase factor associated with the hopping 
is deflned by the gauge potential A^, which we choose to 
be of the following form 



A. 



27rr,- 



(3) 



spanned by all states 



with n'^ + n = h forms the 



where a,j. y are the Pauli matrices, $ is an Abelian flux, 
e-x denotes the unit vector along x, Ti is the spatial co- 
ordinate of site i, and a, j3 are parameters characterizing 
the non-Abelian vector potential. Although interesting 
anisotropy effects can be expected from choosing a ^ /3 
|351 136] ■ in this work we shall consider a — P for sim- 
plicity. With this choice, the intra-species hopping terms 
(aja^ and h\hj) become proportional to cos a, while the 
inter-species hopping terms with non-Abelian vector po- 
tential (a|^bj and h\kj) become proportional to sin a. For 
a = and <I> ^ 0, we thus recover the Hofstadter problem 
of a constant magnetic fleld perpendicular to the two- 
dimensional (2D) system in the Landau gauge [37]. Note 
that also in the opposite limit, <I> = and a 7^ 0, where 
the vector potential becomes constant, the non-Abelian 
character of the gauge potential, i.e. [A^, Aj] 7^ 0, yields 
a constant but finite gauge field. For the Abelian flux, we 
shall focus on $ = p/q, where p and q are co-prime inte- 
gers. Most of our work has been done for $ = or 1/2; 
however, the method developed here can be straightfor- 
wardly extended to other values of $ as shown in Ref. 



III. MOTT INSULATING PHASE 

A. Hopping expansion 

Our approach to study the full Hamiltonian H — 
Hq -|- Hi is based on a resummed hopping expansion as 
developed in Refs. [30] and i31i for the single species BH 
model. In this formalism, one considers the hopping 
term as a perturbation and focuses on the imaginary- 
time evolution of the operators: &{t) 
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Introducing artificial sources jfir), j^ir): Hi [{j}](r) = 
Hi(^) + E^ {jtir)k,{T) + jHtMt) + h.c.) with P the 

inverse temperature, and {j} denoting the set of all four 
currents, the free energy of the system can be written as 
a functional of the sources: 

= -ilnTr(e-^H''T;e-/o''dr H^imKr)^ (4) 

where Tr indicates imaginary-time ordering. From this it 
can be directly seen that the derivatives PSF / 5jl'^^\T) = 
(ai(bi)(T)) — '^'f^\T) yield the order parameter fields 

^^^''■'(t) which vanish within the MI phase. For a de- 
scription of the Mott physics, it is thus sufficient to ex- 
pand up to second order in the currents. Further- 
more, as quantum fluctuations of the hopping term scale 
down with dimension [38^, in a two-dimensional system 
an expansion of up to first order in the hopping 
strength J and a subsequent resummation is expected to 
yield qualitatively reasonable results. This is automati- 
cally achieved by performing a Legendre transformation 
iteratively in the hopping, which substitutes the source 
fields {j} by the physical order-parameter field. Carrying 
out these steps, as detailed in Refs. l30l and [31] we finally 
obtain the effective action of the system up to second 
order in J /U : 



AO 



(5) 



where CJO- = (T,6| (t)6j(t')>o is the unperturbed two- 
point function. Here the operators and Oj may be of 
type a or 6, and the thermal average is with respect to 
Hq. Note that G" is a function of r - r'; thus (G°)~i 
is most easily found in frequency space. We point out 
that an alternative way to deriving the same effective 
action within a random phase approximation is described 
in Ref. E for the BH model, and can straightforwardly 
be generalized to systems with an Abelian gauge field 
[15]. From Eq. we find that G = [{G^r^ - J^J^^ 
is the resummed two-point function which can be used, 
for example, to find the boson momentum distribution 



Uk = — lim — 



(6) 



These momentum distributions can be observed in time- 
of-flight (TOF) measurements, as we shall discuss later. 

The excitation spectra of the bosons can be obtained 
from the poles of G or equivalently by setting up the 
equation of motion: SSH'i'}]^^ /S^'^'^ = 0. In Fourier 
space, this reads 



E 



Sk.k'G'^i'^M) 



Jk.k' 



*fc'(^M) 



0, (7) 



where Jk,k' is the Fourier transform of J^ . After an an- 
alytic continuation to real frequencies iuju uj + ie, the 
solutions of Eq. ([T]) yield the dispersion relations. In the 
following subsections, we apply this general procedure to 
specific choices of gauge potentials and parameters A and 



B. Independent Species 

The simplest non-trivial parameter choice which we 
shall treat in this section corresponds to A = and 
< n/U < 1. In this case, Hq describes two independent 
standard BH systems with unique non-degenerate ground 
state having one boson of each species per site. The 
coupling between them is provided by the inter-species 
hopping terms arising from the non-Abelian gauge po- 
tential. For this case, it is clear that the unperturbed 
Green function has vanishing off-diagonal terms, i.e. 
G?2 - (a^(0)b(T))o = and - (bt(0)a(T))o = 0. 
From our symmetric choice of parameters, it also fol- 
lows that G?i = G^2 rendering G° 1 2x2- Further, the 
site-factorizable nature of Hq guarantees that the diag- 
onal elements of G'^ are given by 6k,k' times a function 
G^{ujm) of a single Matsubara frequency 



G°(c^m) = 



n + 1 



n,m— 



(8) 



where Z° = Li^^m^o e''^"'" ' ^» = -^"."i - En-\^m- 
In the zero-temperature limit, the Boltzmann sums in Eq. 
([S]) reduce to a single term corresponding to the ground- 
state occupation numbers, n = m = \. In the following 
two subsections, we shall compute G for the two simplest 
choices of gauge potentials corresponding to a constant 
non-Abelian gauge field: the one without Abelian flux, 
$ = 0, and that with an Abchan flux $ = 1/2. 



1. $ = 

For a constant gauge potential without Abelian flux, 
$ = 0, the hopping matrix is diagonal in momentum 
space, i.e., 

JkM' =2J{ cosa[cos(fca;) -|- cos(fcy)]l (9) 
- sina[sin(A:^)aj^ + sm{ky)a^\\^k,k' ^ 

where here and in the rest of the paper, we have set the 
lattice spacing a = 1. As the unperturbed Green function 
G° is already diagonal, we need to diagonalize only the 
hopping matrix. This yields 



cos(a) [cos(fc:i;) -t- cos(fcj,)] 
± sin(a)y^cos2(fc:r) 4-cos2(fcy) 



(10) 
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FIG. 1. Boson energy dispersion for ^ — 0, a = 1, jj,/U 
0.41, and J/U = 0.05. See text for details. 



The energy bands of the system at zero temperature are 



thus given by [G"(w + iry)]"i|T=o - 
have Wick rotated back to real frequency, 
two quadratic equations 



- uj{2ii 

- u;(2/i - 



U 

U ^ 



0, where we 
This yields 







(11) 



leading to four energy bands shown in Fig. [T] Two of 
these bands occur at w > and the other two at a; < 0, 
so that they may be interpreted as particle/hole excita- 
tion spectra of the system. We note that as a — > 0, the 
particle and the hole spectra become increasingly similar 
and ultimately indistinguishable for a — Q yielding the 
standard dispersion of MI bosons with no gauge poten- 
tial 0. Furthermore, the particle and hole excitations 
have a gap in the Mott phase which closes as one ap- 
proaches the superfluid-insulator transition by increas- 
ing J /U . Beyond the transition point, which occurs at 
J = Jc, the solutions of Eq. ( 11 ) are complex, indicating 
instability of the Mott phase. 

The superfluid-insulator phase boundary, as obtained 
by the procedure described above, is shown in Fig. [2] We 
find that the phase boundary has the usual lobe struc- 
ture. However, the value of Jc at the tip of the Mott 
lobe is strongly influenced by a and can thus be tuned 
by varying the strength of the gauge field. This leads 
to re-entrant superfluid-insulator transitions in the sys- 
tem by variation of a, provided that J is appropriately 
fixed at (say) J = Jc(a — 0.2tt) as can be seen from 
Fig. [2] We also note that the universality class of the 
superfluid-insulator transition has the same properties 
as in the standard BH case [1]. At the lobe tip, the 
additional particle-hole symmetry renders the dynamical 
critical exponent z of the transition to be unity; at other 
points, z — 2. 

One of the key difference of the superfluid-insulator 
transition in the present system from the normal BH 
model is that the position of the minima of the low-energy 
excitations at and near the critical point in the first Bril- 
louin zone are at finite momenta and strongly depend on 
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FIG. 2. Up The Mott lobes with a dependent heights. The 
solid (dashed) lines indicate the lobes for a = 0(1). Down: 
Value of Jc at the tip of the n = 1 Mott lobe as a function of 
a. The sohd (dashed) lines correspond to $ = 0(1/2). 




a/ n 



FIG. 3. The position of the boson energy minima at finite 
wavevector parametrized by k (see text) in systems without 
flux (solid line) and with 1/2 flux (dashed line). In the latter, 
a QPT of the excited band occurs at a = n/4. 



a (see Fig. [s]). Since the position of these minima cor- 
respond to precursor peaks of the bosons near the crit- 
ical point [7], this feature is reflected in the momentum 
distribution of the bosons in the Mott phase near the 
transition. From the Green function G, we compute the 
momentum distribution using Eq. ^. The results of this 
calculation are shown in the top panel of Fig. |4j We find 
that at finite a, the precursor peaks of the momentum 
distribution at J ~ 0.97 Jc{a) are at finite momenta, re- 
fiecting the fact that the subsequent condensation of the 
bosons aX J — Jc{a) will occur at non-zero momenta. We 
note that such a pattern in the momentum distribution 
should be easily picked up in TOE experiments. Since 
standard TOE experiments will measure the distribution 
of both species simultaneously, the presence of the non- 
condensing mode will slightly reduce the visibility of the 
pattern. However, we expect that one should be able to 
easily subtract this background in order to observe the 
sharp peaks of the condensing mode. The most signifi- 
cant feature distinguishing the non-Abelian scenario from 
the known standard one, is the number of peaks which 
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in general is quadrupled by the non-Abelian gauge field. 
We predict this feature to be clearly observable in TOF 
experiments. 



2. $ = 1/2 

Turning to the case with an Abelian flux $ in the gauge 
potential of Eq. ([s]) , we find the general structure for the 
hopping matrix to be given by: 



(12) 



making Eq. ([t]) to have off-diagonal terms connecting 
momenta in the magnetic Brillouin zone which differ by 
±27r<I>. Here A, B and C are functions of momentum 
which will be specified later. For a generic $ = p/q, the 
periodicity of the lattice '^k = '^k+2T!ne^ ensures that Eq. 
([t]) leads to a set of q closed equations. To write them 
down, we introduce the notation ^fc+„x27r$e^ — ^k.n 
with n — 0, ■ ■ ■ , q — 1. In this notation, the equation of 
motion can be written in the Harper-like form: 



with 

M(fc„n) = [G"(c^m)]"'12x2-2J 

/ cos(afc2: -I- 27r<i>n) cos a ism{akx 
\ — i sin(afcx + 27r$n) sina cos(afc3; 



-1=0, 
(13) 



(14) 



- 27r<i>n) sin a 
27r<i>n) cos a 



and 



TV EE -J 



cos /3 i sin /3 
i sin (3 cos (3 



(15) 



As each equation in Eq. ( |13[ ) is a two-component equa- 
tion, for a flux of $ = p/q, we have 2q closed equations. 
Here we shall focus on p = l,q = 2 which allows us to 
find the solutions of these equations analytically via di- 
agonalization of a 4 x 4 hopping matrix. 

We find that the presence of the magnetic fiux $ = 1/2 
splits each band, so that we now have four particle and 
four hole excitations. Again the most relevant bands 
are the particle (hole) excitation at lowest (highest) fre- 
quency. The gap between these excitations closes at 
J = Jc leading to a second-order QPT separating the 
MI and the SF phases. The lobe structure of the phase 
boundary (and also the universality class of the transi- 
tion) remain unchanged by the flux, as can be seen from 
Fig. [2] We note however that the plot of Jc(a) as a func- 
tion of a, shown in the bottom panel of Fig. [2j has a qual- 
itatively different behavior compared to the case with- 
out Abelian flux discussed previously; nonetheless, the 
system will show similar re-entrant superfluid-insulator 
transitions as a is varied for a fixed J. 
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FIG. 4. Momentum distributions of the bosons at J « 
0. 97 Jc. Bright regions correspond to high densities. The 
upper (lower) rows correspond to $ = 0(1/2). The non- 
Abelian field strengths a are, from left to right in each panel, 
0, 0.7, 0.8, and7r/2. 



The most interesting difference between $ = and 
$ = 1/2 concerns the number and positions of the ex- 
trema of the particle and hole excitations. We find 
that within the first magnetic Brillouin zone (k^ € 
[— 7r/(7, 7r/q], ky G [— 7r,7r]) and in the absence of a 
non-Abelian field (a = 0), there are two extrema at 
k — (0, 0), (0, ±7r), in agreement with Ref. 25 Denoting 
these three points in the Brillouin zone as fco; 's+i ^^nd 
fc_ , the extrema for non-zero a can be shown to occur at 
feo + (±fc, ifc), ^+ + (ifc; —k), and fc_ — {±k, +k), where 
fc as a function of a is plotted in Fig. [3j From this plot, 
we find that as long as a < 7r/4, we have eight extrema. 
When a = 7r/4, we get fc = 7r/2 and the extrema are com- 
pletely shifted to the zone edges, i.e. again we have only 
two extrema per Brillouin zone. For 7r/4 < a < 37r/4, 
a plateau with a single fixed extremum is found. The 
derivative dk/da diverges at a — >■ 7r/4 which means that 
at this critical value of the non-Abelian field the global 
minima of the excited bands abruptly change their posi- 
tion. Thus a slight change in the gauge field strength is 
expected to completely modify the dynamical behavior 
of the system. This remarkable behavior can be directly 
observed in the momentum distribution shown in the bot- 
tom panel of Fig. [4j In particular, the abrupt change in 
the pattern of the momentum distribution when the flux 
is varied across 7r/4 reflects the sudden change in the po- 
sition of the band minima with small change in a. This 
behavior is reminiscent of the QPT of excited states dis- 
cussed in Ref. [311 Also the abrupt sign reversal of the 
Hall conductivity at half filling in a system of hard-core 
bosons subjected to a gauge field as studied in Ref. [33| 
falls into this category of phenomena where some control 
parameter modifies the system's dynamics in a discon- 
tinuous way. 




FIG. 5. The excitation spectrum of the bosons in the XY 
phase system for a = 1, $ = 0, ^/U = 0.15, A = 0.3, and 
J/U = 0.021. 



C. XY configuration 

In this section, we consider the case where A > such 
that the ground state of the two-species model without 
the gauge field correspond to the XY phase discussed 
in Ref. [121 where the bosons are in the superposition 
of both the states. Consequently, will have non-zero 
off-diagonal components and one needs to diagonalize the 
full matrix (G°)~^ — Jfc in the presence of the non-Abelian 
flux. 

We begin by computing the elements of (7° at T = 
which are given by 



AT.aI(r)a,(T'))GS (T.aI(r)b,(r'))GS 
^(T.bl(T)a,(r'))GS (T.bl(r)b,(r'))GS 



5t/ 



(16) 



For the XY ground state |GS) = ^(|10) |01)), dis- 
cussed in ReffT^ this matrix reads 



1 



1 



ioj) [/ — /i - 
UX 



2{n + iuj){-U\ + ^ + iw) 



2iuj 
(17) 



1 



where 1 is the unit matrix and ax denotes the Pauli ma- 
trix. Inserting Eq. (17) into Eq. ([7]), and considering the 
case $ = 0, we find that the presence of the off-diagonal 
elements in leads to two independent equations for 
the band dispersions given by 



Ml + Ak + jAq + iBfeP + AhiBk + Bl)=Q 



Ml +Ak^ JNq + + M2{Bk + Bl) = (18) 



FIG. 6. Left: The critical hopping Jc marking the tip of the 
lobe for A = 0.3 (solid line) and A = 0.6 (dashed line) as a 
function of a. Right: The minima (maxima) of the lowest 
(highest) particle (hole) excitation are found at kya = —a 
and k = ±kx as shown in the plot (the solid line corresponds 
to A = 0.3, the dashed line to A = 0.6). It is different from 
or ±7r only within a small region around a — it/2, in which 
the single extremum suddenly splits into two. 



where 



Ml = 



M, = 



{Lu + fi- U){U^X + UX{uj + 2{n + uf) 

2C/2A - 2(/x w)2 
UX{^i + uj- Uf 



2U^X-2{n + uj)^ 
Ak = 2 J cos a(cos + cos ky) 
Bk = 2Jsina(sin/cj, -I- isinfe^;) 



(19) 



Solving the first of these equations, we find two pos- 
itive {E > 0) and one negative {E < 0) solutions, 
while the second equation has solely one positive solution 
which is the second highest band. The resulting band- 
structure is shown in Fig. [Sj Note that here the pres- 
ence of the off-diagonal component of the Green function 
which originates from the XY ground state leads to more 
particle-like than hole-like excitations. This feature is a 
consequence of particle-hole asymmetry originating from 
jjah _^ jj ^ Also, the splitting of the two highest particle- 
like excitations is a consequence of the non-Abelian na- 
ture of the hopping. This splitting vanishes in the limit 
a ^ 0. 

For the QPT into the SF state these higher modes 
again do not play a role. The Mott lobe, on which at 
least one mode becomes gapless, now extends from < 
jjL < XU . The value of Jc marking the height of the 
lobe depends on both a and A, as illustrated on the left 
panel of Fig. |6] Furthermore, as in the cases studied 
before, we find that the minima of the dispersion occur 
at finite wavevectors. However, in contrast to the cases 
studied before, the band spectrum in the XY phase is 
not symmetric under ky — )> 



-ky. This property of the 



dispersion can be traced back to G since Bk in Eq. ( 19 ) 
is not invariant under such a transformation. 

The above-mentioned properties of the dispersion of 
the particle- and hole-like excitations are reflected in the 
nature of the momentum distribution of the bosons in 
the Mott phase near the quantum critical point. We find 
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1-1 

kya/n 



1-10 1 



FIG. 7. Momentum distributions at J « 0. 97 Jc in XY 
system for A — 0.3. The bright regions correspond to high 
densities. The field strengths a are, from left to right, 
0, 1, 1.4, 1.5, 7r/2, 1.65, 1.75, andyr. 




1-1 1-1 1-1 

kya/n 



FIG. 8. The momentum distributions of the bosons at J ~ 
0.97Jc in the XY system with ^ = 1/2 for A = 0.3 and field 
strengths a, from left to right, 0, 0.7, 1, and 1.5. 



that for any given a, the condensing modes are located 
at ky = —a. As shown in Fig. [Tj there are no analogous 
peaks at ky = a which clearly reflects the breaking of 
the parity symmetry discussed above. Also we note that 
for the range < a < 1.5 (1.65 < a < tt), there is a 
single condensing mode at kx — [k^ = ±7r). Around 
a = 7r/2, however, the condensing mode splits into two 
at momenta k — (±fc, — a), where k is plotted on the 
right panel of Fig. |6] as a function of a. We note that 
such a splitting can modify the dynamical behavior of 
the system. 

Finally, we study the influence of a magnetic flux on 
the XY system, as done before for A = 0. As expected 
we find the Abelian flux $ = 1/2 to split each band into 
two, such that the system exhibits two hole and six par- 
ticle excitations. As before, for a = 0, the presence of 
this flux provides the band structure with two extrema 
at fe = (0,0) and (0, ±7r). The resulting momentum dis- 
tribution near condensation is shown in Fig. [Sj As a is 
increased, the position of the peaks of the momentum 
distribution initially shifts along ky. However, close to 
a = 7r/2, the peaks split along kx leading to four peaks 
which finally reach the zone edge k — (±7r/2, ±7r/2) at 
a = tt/2. These features, as in the case of $ = 0, should 
be experimentally observable via standard TOF experi- 
ments. 



IV. SUPERFLUID PHASE 



that, as pointed out in Ref. [7] and [551 the quadratic ac- 
tion (Eq. ([5])) does not capture the physics of the ordered 
phase but needs to be supplemented by the quartic term 
at the mean-field level. These terms can be directly com- 
puted using the methods of Ref. [71 and [251 within a strong 
coupling expansion; however, it is often enough to guess 
their forms from the symmetry of the underlying Hamil- 
tonian. In this section, we are going to take the latter 
route and chart out the characteristics of the resulting 
SF phase. 

For determining the order parameters in the SF phase 
we need to construct the quartic part of the effective ac- 
tion in terms of the order parameter fields and minimize 
it. To this end, we first rewrite the quadratic action by 
diagonalizing its kernel as shown in Section |III[ Let us 
denote the eigenvalues and eigenmodes of the quadratic 
action as mn^k{^) and 5'„^fc((u;), respectively. In the case 
when the Abelian fiux is $ = p/q, n varies from to 
2(7—1. Note that here we have adopted the convention 
that for $ = 0, the system remains with only two eigen- 
modes, so that we have n = 0, 1. In the zero-temperature 
limit we may then write, 

2g-i „ 

•^^^ = E E / ^nM\^n.k{^)\'- (20) 

Ti=0 k •' 

Here, the sum over k is restricted to the first magnetic 
Brillouin zone. The fourth-order term can be written in 
this basis as 

29-1 i-f) 

- .9/2 Y.Y.h^^ |§„(r„T).M'„(n,T)p(21) 

ri=0 i ''^ 

where we have transformed the '^n and fields to real 
space. Here, g > is the exact two-particle vertex func- 
tion of the bosons in the local limit, which has been com- 
puted in Refs. [3 |3T] and i denotes lattice sites. 
With this the SF action may be written as 

^SF^^MI^^(4)^ (22) 

Now we note that at the onset of superfiuidity only one 
of the eigenmodes condenses; therefore, it is possible to 
analyze the SF phase within a mean-field approximation 
by ignoring the other modes. This observation allows us 
to get rid of the n-sum over all bands in Eq. (21): 

5(4)^5/2^ /%T |i'(r„r).vl,(r„r)|2. (23) 
I Jo 

In the presence of a non- Abelian gauge field with mag- 
netic fiux <I> = 0, or flux $ — p/q = 1/2, we have one 
or more than one minima of the boson energy spectrum 
depending upon the non- Abelian fleld strength a. If the 
particle/hole modes consist of s degenerate minima, then 
the corresponding Ginzburg-Landau theory can be ex- 
pressed by s low-energy fluctuating two-component fields 
(order parameters) (r, t) around these minima [251 126) : 



In this section, we investigate the nature of the SF 
phase into which the transition takes place. We note 



*(n,t) = ^x„(ri),0„(ri,t). 



(24) 



where we have Wick rotated to real time. The coefficients 
xifi) are the real space eigenfunctions corresponding to 
the minimum energy band at fc = A:™™) which can 

be expressed as: 



x«(n) 
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5]c,e(< 

1=0 



(25) 



Note that the sum in the above expression is restricted 
to q terms, since the functions Xnifi) describe only the 
part of the spatial dependence of \E'(ri,t) that can be 
factored out for each term in the sum in Eq. ( |24| ). Here 
ci denotes the components of eigenvectors corresponding 
to the minimum energy band at fe = (fc™'", fc™™)- 

In general, the quartic part of the Landau-Ginzburg 
action is difficult to obtain, since it is restricted only by 
the invariance under projective symmetry group (PSG) 
of the underlying square lattice [3S]. The elements of 
PSG include in this case translations along the x and y 
axes, rotation by tt/2 around the z axis, and reflections 
about X and y axes. In our case, the situation is, however, 
much simpler, since we know the microscopic form of the 
quartic action, Eq. (21). We may therefore substitute 



Eq. (24) into Eq. (21 ), and obtain the explicit form of the 



quartic action in terms of the order parameters (r^ , t) . 
We can then find the saddle point of the total action with 
respect to (f>n {vi , t) , and thus directly calculate the order 
parameters in the SF phase. 

Let us first consider the case where the number of both 



flavors at each site is equal, as discussed in Sec. Ill B If 



the boson spectrum has one minimum in the magnetic 
Brillouin zone, then the corresponding low-energy fleld 
can be written as 



'^ir^,t) = xi('^i)0i(n,t), 
and thus the SF density reads 



9-1 
1=0 



(26) 



(27) 



For $ = 0, we can get rid of the sum in Eq. ( 25 1 and ps is 



simply equal to \4>i\ , which has no modulation along x. 
In contrast, for $ = p/q = 1/2, we find that ps displays 
a spatial pattern. 

Next, we consider the case where there are two min- 
ima at kf" = {n/2, n/2) and fcf'™ = (7r/2, ~n/2) within 
the magnetic Brillouin zone for a — tt/A and $ = 1/2. 
Note that these are on the zone edge, such that the min- 
ima at the opposite edge are equivalent. In this case, 
'0(n,t) = Xllrl)(|)l{r^,t) + X2{ri)(|)2{r^,t). Following 
the coarse-graining procedure charted out in Ref. we 
find that the SF ground state corresponds to the con- 
densation of any one of the low-energy fluctuating fields 
^ 0, (02) = or ((/.i) = 0, (02) ^ 0. The corre- 
sponding plot for a = 7r/4 and $ = 1/2 in Fig. [9^, shows 
a similar periodic pattern as found for a single minimum 
for $ = 1/2. Similar analysis can be done for four min- 
ima at (±7r/2, ±7r/2) for a = tt/2 and $ = 0; in this case 





FIG. 9. a) Plot of SF density p^/po for <!> = 1/2 and a = 7r/4 
at fe""'" = (7r/2,±7r/2) for even filling, b) Same plot for $ = 
1/2 and a = tt/2 at fc""" = {tt/2, ±tt/2) for XY phase. 



we find that only one out the four field condenses; conse- 
quently there is no modulation of SF density. Note that 
we have restricted ourselves so far to the minima occur- 
ring at the wavevectors (tt/si, tt / S2), where si^2 = il- In 
principle the analysis can be extended to the situations 
when the minima occur at (777/2, 5tt /2) with rational and 
small 7 and 5] however, the analysis becomes technically 
involved and we have not attempted that in this work. 
We do not have general understanding of implementing 
the above-mentioned coarse-graining procedure for irra- 
tional 7, 5. 

Finally, we briefly comment on the SF density in the 
XY phase. Following the procedure discussed before, we 
again find a constant SF density for any non-Abelian 
gauge field with $ = 0. For a = tt/2 and $ = 1/2, there 
are two minima of the spectrum, and we find that the 
SF ground state corresponds to the condensation of both 
fields around these minima. The corresponding plot for 
a — TT /2 is shown in Fig. [Oja. 

Thus we generically find that in the presence of a non- 
Abelian gauge field, the SF density displays a spatial pe- 
riodic pattern if there is a finite flux (<& — 1/2); however, 
there is no such modulation without flux, $ = 0. The 
method that has been discussed here can be used for any 
filling fraction ^ = p/q. We expect different spatial pat- 
terns of the SF density for other $, and leave the detailed 
analysis of it for future study. 
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V. CONCLUSION 



In conclusion, we have studied the Mott phases and 
the superfluid-insulator transition of strongly-interacting 
two-species bosons in the presence of a non-Abelian 
gauge field. We have shown that such a system of bosons 
is expected to display novel patterns of precursor peaks 
in the momentum distribution in the Mott phase close 
to the superfluid-insulator critical point. We have also 
demonstrated the presence of a re-entrant superfluid- 
insulator transition as the strength of the non-Abelian 
field is varied for a fixed hopping amplitude. Finally, we 
have found that the presence of an additional Abelian 
field with $ = 1/2, leads to spatial modulation of the su- 
perfluid density in the SF phase near the critical point; 
however, no such modulation is expected for $ = 0. 
We note that the precursor peaks in the momentum dis- 



tribution and the presence of the re-entrant superfluid- 
insulator transition can be easily detected in standard 
TOF experiments. Finally, we point out that we have 
restricted ourselves to model parameters for which the 
Mott phase has translational symmetry; it will be inter- 
esting to extend our study to the case where the Mott 
state has broken translational symmetry. We leave this 
as a subject of future study. 
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